Complementary descriptions of shape/phase 
transitions in atomic nuclei 

Dennis Bonatsosi], E. A. McCutchanH, N. V. Zamfiril 

" Institute of Nuclear Physics, N.C.S.R. "Demokritos" , 

GR-15310 Aghia Paraskevi, Attiki, Greece 
^ Wright Nuclear Structure Laboratory, Yale University, 
New Haven, Connecticut 06520-8124, USA 
QQ ", '^ National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania 

O 

o 

^ ' Abstract 



00 






Shape/phase transitions in atomic nuclei have first been discovered in the framework of 
the Interacting Boson Approximation (IBA) model. Critical point symmetries appropriate 
CS ! for nuclei at the transition points have been introduced as special solutions of the Bohr 
Hamiltonian, stirring the introduction of additional new solutions describing wide ranges 
of nuclei. The complementarity of the IBA and geometrical approaches will be demon- 
strated by three examples. First, it will be shown that specific special solutions of the 

O I Bohr Hamiltonian correspond to the borders of the critical region of the IBA. Second, it 

will be demonstrated that the distinct patterns exhibited in different transitional regions 
by the experimental energy staggering in 7-bands can be reproduced both by the IBA and 
by special solutions of the Bohr Hamiltonian. Third, a first attempt to obtain a IBA SU(3) 

>• ' level scheme from a special solution of the Bohr Hamiltonian will be presented. 

O 

^ : 1 Introduction 

O \ Atomic nuclei are known to exhibit changes of their energy levels and electromagnetic tran- 

sition rates among them when the number of protons and/or neutrons is modified, resulting 
in shape phase transitions from one kind of collective behaviour to another. These transi- 
tions are not phase transitions of the usual thermodynamic type. They are quantum phase 

/\ \ transitions [1] (initially called ground state phase transitions [2l|), occurring in Hamiltonians 
c^ ■ of the type H = c{Hi + gH2), where c is a scale factor, g is the control parameter, and Hi, 
H2 describe two different phases of the system. The expectation value of a suitably chosen 
operator, characterizing the state of the system, is used as the order parameter. 

In the framework of the Interacting Boson Approximation (IBA) model [3], which de- 
scribes nuclear structure of even-even nuclei within the U(6) symmetry, possessing the 
U(5), SU(3), and 0(6) limiting dynamical symmetries, appropriate for vibrational, axially 
deformed, and 7-unstable nuclei respectively, shape phase transitions have been studied 25 
years ago [2] using the classical limit of the model [H El [6], [7] , pointing out that there is (in 
the usual Ehrenfest classification) a second order shape phase transition between U(5) and 
0(6), a first order shape phase transition between U(5) and SU(3), and no shape phase 
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transition between 0(6) and SU(3). It is instructive to place |T] these shape phase tran- 
sitions on the symmetry triangle of the IBM |8], at the three corners of which the three 
limiting symmetries of the IBM appear. 

More recently it has been realized [01 [ID] that the properties of nuclei lying at the 
critical point of a shape phase transition can be described by appropriate special solutions 
of the Bohr Hamiltonian [IT], labelled as critical point symmetries. The E(5) critical point 
symmetry [9] has been found to correspond to the second order critical point between U(5) 
and 0(6), while the X(5) critical point symmetry [TU] has been found to correspond to the 
first order transition between U(5) and SU(3). 

The introduction of the critical point symmetries E(5) [9j and X(5) [10] has triggered 
much work on special solutions of the Bohr Hamiltonian, corresponding to different physical 
situations [121 [13] • Several of these solutions will be mentioned below, where the comple- 
mentarity of the IBA and geometrical approaches will be demonstrated by three examples, 
given in Sections 2, 3, and 4. 

2 Connecting the X{5)-f3\ X(5)-/3^ and X(3) models 
to the shape/phase transition region of the Inter- 
acting Boson Model 

As already mentioned in the introduction, shape/phase transitions in atomic nuclei were 
first investigated [2J within the Interacting Boson Approximation (IBA) model \3\ by con- 
structing the classical limit of the model, using the coherent state formalism [H [5l [6l [7] . 
Using this method it was shown [21 [B] that the shape/phase transition between the U(5) 
[spherical, vibrational] and SU(3) [prolate axially symmetric deformed, rotational] limiting 
symmetries is of first order, while the transition between the U(5) and 0(6) (7-unstable) 
limiting symmetries is of second order. Furthermore, the region of phase coexistence within 
the symmetry triangle [8] of the IBA has been studied [TU [T5| [T6] and its borders have been 
determined [171 [18] , while a similar structural triangle for the geometric collective model 
has been constructed [T9] . 

Recently, special solutions of the Bohr Hamiltonian, called critical point symmetries 
[9l [To] , describing nuclei at the points of shape/phase transitions between different limiting 
symmetries, have attracted considerable attention, since they lead to parameter indepen- 
dent (up to overall scale factors) predictions which are found to be in good agreement 
with experiment [201 [211 [221 123]- The X(5) critical point symmetry [TU], was developed 
to describe analytically the structure of nuclei at the critical point of the transition from 
vibrational [U(5)] to prolate axially symmetric [SU(3)] shapes. The solution involves a 
five-dimensional infinite square well potential in the (3 collective variable and a harmonic 
oscillator potential in the 7 variable. 

The success of the X(5) model in describing the properties of some nuclei with parameter 
free (except for scale) predictions has led to considerable interest in such simple models to 
describe transitional nuclei. Since its development, numerous extensions involving either no 
free parameters or a single free parameter have been developed. Those approaches which 
involve a single parameter include replacing the infinite square well potential with a sloped 
well potential [21] , exact decoupling of the /? and 7 degrees of freedom [25] , and displacement 



of the infinite square well potential, or the confined /3-soft model ^26j. Parameter free 
variants of the X(5) model include the X(5)-/?^ and X(5)-/9^ models [27j, in which the 
infinite square well potential is replaced by a (3"^ and a /5^ potential respectively, as well as 
the X(3) model [28], in which the 7 degree of freedom is frozen to 7 = 0, resulting in a 
three-dimensional Hamiltonian, in which an infinite square well potential in (3 is used. 

It is certainly of interest to examine the extent to which the parameter free (up to overall 
scale factors) predictions of the various critical point symmetries and related models, built 
within the geometric collective model, are related to the shape/phase transition region 
of the IBA. It has already been found [29] that the X(5) predictions cannot be exactly 
reproduced by any point in the two-parameter space of the IBA, while best agreement is 
obtained for parameters corresponding to a point close to, but outside the shape/phase 
transition region of the IBA. 

In a recent paper pO] the parameter independent (up to overal scale factors) predictions 
of the X(5)-/3^, X(5)-/9^, and X(3) models, which are variants of the X(5) critical point 
symmetry developed within the framework of the geometric collective model, are compared 
to the results of two-parameter Interacting Boson Approximation (IBA) model calculations, 
with the aim of establishing a connection between these two approaches. It turns out 
that both X(3) and X(5)-/5^ lie close to the U(5)-SU(3) leg of the IBA symmetry triangle 
and within the narrow shape/phase transition region of the IBA. In particular, X(3) lies 
close to Ccrit [HI 5 the left border of the shaded shape/phase transition region of the IBA, 
corresponding to IBA total energy curves with two equal minima, while X(5)-/9^ lies near 
the right border of the shape/phase transition region, C,** [17] , corresponding to IBA total 
energy curves with a single deformed minimum. A set of neighboring even-even nuclei 
exhibiting the X(3), X(5)-/5^, and X(5)-/3^ behaviors have been identified (^^^Os-^^^Os- 
i^'^Os) I30j. Additional examples for X(3), X(5)-/32, and X(5)-/5^ are found in ^^^Pt, ^^^Ce, 
and ^^^Er, respectively [30]. The level of agreement of these parameter free, geometrical 
models with these candidate nuclei is found to be similar to the predictions of the two- 
parameter IBA calculations. 

It is intriguing that the X(3) model, which corresponds to an exactly separable 7-rigid 
(with 7 = 0) solution of the Bohr collective Hamiltonian, is found to be related to the IBA 
results at Ccrit, which corresponds to the critical case of two degenerate minima in the IBA 
total energy curve [18], approximated by an infinite square well potential in the model. It 
is also remarkable that the X(5)-/3^ model, which uses the same approximate separation of 
variables as the X(5) critical point symmetry, is found to correspond to the right border 
{(,**) of the shape/phase transition region [17], related to the onset of total energy curves 
with a single deformed minimum, comparable in shape with the (3"^ potential used in the 
model in the presence of a L{L + l)/(3/3^) centrifugal term [25] . 

Comparisons in the same spirit of the parameter independent predictions of the E(5) 
critical point symmetry [9] and related E(5)-/5^" models [ST, ^32], as well as of the related to 
triaxial shapes Z(5) [33] and Z(4) ^34j models, to IBA calculations and possible placement 
of these models on the IBA-1 symmetry triangle (or the IBA-2 phase diagram polyerdon 
[35l [361 EIj) can be illuminating and should be pursued. 

It should be noticed that Ref. [30] has been focused on boson numbers equal or close 
to 10, to which many nuclei correspond. A different but interesting question is to examine 
if there is any connection between the X(3), X(5)-/3^, and/or X(5)-/5'^ models and the 
IBA for large boson numbers. This is particularly interesting especially since it has been 



established (initially for N = 1,000 [Mj, recently corroborated for N = 10,000 ^) that 
the IBA critical point of the U(5)-0(6) transition for large A^ corresponds to the E(5)-/3^ 
model, i.e. to the E(5) model employing a /?^ potential in the place of the infinite well 
potential 



3 Staggering in 7 bands and the transition between 
different symmetries of nuclear structure 

The importance of the 7 degree of freedom in nuclei with static quadrupole deformation 
has been known for decades starting with the work of Wilets and Jean [39j , and Davydov, 
Filippov, and Chaban ^01 HI] . Since the bandhead of the quasi-7 band is not fixed in most 
geometrical models, we investigate instead, the parameter free predictions of the spacings 
within the quasi-7 band. For this purpose, we use the odd-even staggering in gamma bands 
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[E{J^) - E{{J - 1)+)] - [E{{J - 1)+) - E{{J - 2)+) 



E{2t) 

which measures the displacement of the {J—l)'t level relative to the average of its neighbors, 
J^ and (J — 2)j;, normalized to the energy of the first excited state of the ground state 
band, 2^. Since S{J) is a (discrete) derivative, it is very sensitive to structural changes. For 
example, the energy levels of the 7 band in a 7-independent potential [35] cluster as (2+), 
(3+,4+), ..., opposite to the rigid triaxial rotor (2+, 3+),(4+, 5+), ...clustering pattern 



In a recent paper |33] , three categories of transitional regions have been considered: 

1) The 7-soft region between the vibrator and a deformed 7-soft structure where the 
potential is 7-independent. This corresponds to the U(5) to 0(6) transition in the language 
of the Interacting Boson Approximation (IBA) [3]. This is the region containing the critical 
point symmetry E(5) [9], as well as the second-order phase transition in the IBA between 
U(5) and 0(6) [21 [6]. The Xe, Ba, and Ce series of isotopes provide good manifestations 
of this region [21 US]- It is seen that S{J) exhibits strong staggering with minima at even 
J and maxima at odd J . This behavior is reproduced in the geometrical framework by the 
chain of parameter-independent models formed by E(5) [H], in which an infinite well u{(3) 
potential is used, and the E(5)-/32" {n = 1, 2, 3, 4) models [3ll[32], in which m(/5) = /52"/2. In 
all cases the potential is 7-independent. In the IBA it corresponds to x = and increasing 

c[i5i[n]. 

2) The axially 7-rigid region between the vibrator and the axially symmetric rotor, 
characterized by a harmonic oscillator in 7 with the minimum in 7 close to zero. This is 
the U(5) to SU(3) transition region of the IBA |3], in which a first-order phase transition 
occurs [21 E]- This is also the region where the critical point symmetry X(5) [TD] is found. 
Several Sm, Gd, Dy, Er, U, and Fm isotopes exhibit this behavior ^ |13] . It is seen that 
S{J) exhibits weak staggering with minima at even J and maxima at odd J. This behavior 
is reproduced in the geometrical framework by the chain of parameter-independent models 
formed by X(5) [10], in which an infinite well u{f3) potential is used, and the X(5)-/?^"' 
(n = 1, 2, 3, 4) models [2Z1, in which u{P) = /?^"'/2, with a potential of the form m(/3) +f (7) 



used in all cases, as well as by the exactly separable ES-X(5) and ES-X(5)-/3^ models [H], in 
which potentials of the form m(/5) + f (7)//?^ are used. In all cases f (7) is a steep harmonic 
oscillator potential centered at 7 = 0°. In the IBA this region corresponds to x = —1-32 
and increasing ( [151 IH] • 

3) The triaxial 7-rigid region between the vibrator and the rigid triaxial rotator [iQlHT], 
characterized by fixed 7- values between 0° and 30°, which has no analog in the framework 
of IBA-1 [3]. Few nuclei (^^^Ru, ^™Et, ^^^Os, ^^^Pt, ^^^Th) are known to show this behavior 
[13 J. in which S{J) exhibits strong staggering with minima at odd J and maxima at even 
J. This behavior is reproduced in the geometrical framework by the chain of parameter- 
independent models formed by Z(5) [33J, Z(5)-/3^, Z(4) [34], and Z(4)-/3^, which are special 
solutions of the Bohr Hamiltonian for 7 = 30° in five and four dimensions respectively. 

Having discussed the dependence of S'( J) on angular momentum, it is interesting to focus 
attention on S'(4), the relative displacement of the 3^ state relative to the average of its 
neighbors, 2j; and 4+ normalized to the energy of the 2^ state, and discuss its variation as a 
function of changing structure, using as a structure indicator the ratio -R4/2 = E{4:i)/E{2i). 
The following results are found [43]: 

i) In the U(5)-0(6) region the IBA predicts >S'(4) decreasing with increasing -R4/2. The 
same prediction is made by the geometrical models mentioned in 1), completed with 0(5)- 
CBS [15], which corresponds to E(5) with an infinite well potential displaced from the 
origin. The Xe and Ba chains of isotopes provide an experimental manifestation of this 
behavior. 

ii) In the 0(6)-SU(3) region the IBA predicts S'(4) increasing with increasing -R4/2. The 
same prediction is made by the geometrical models mentioned in 3). 

iii) In the U(5)-SU(3) region the IBA predicts S'(4) decreasing with increasing _R4/2 
for low values of -R4/2, then abruptly jumping to a value close to zero, close to which it 
remains with further increase of -R4/2. The geometrical models mentioned in 2) reproduce 
the same behavior (slight increase of S'(4) with increasing -R4/2) beyond the abrupt change, 
a situation corroborated by experimental evidence from the Nd, Sm, Gd, Dy, and Er isotope 
chains [13]. The abrupt change appears to be related to the first order phase transition 
[21 [6] occuring in this region. It appears that phase transitional behavior occurs close to 
where 5'(4) crosses zero. 

4 Occurence of IBA SU(3) degeneracy in the geomet- 
rical model 

A long standing problem is the derivation from the Bohr Hamiltonian \TT] of a spectrum 
similar to that of the SU(3) limit of the Interacting Boson Approximation (IBA) model [3]. 
The main features of the spectrum should be: 

a) The energy spacings among the 2"^, 4"*", 6"^, ...levels within the ground, (3 and 7 
bands should be identical. 

b) Furthermore, the 2+, 4+, 6+, ... levels of the /? and 7 bands should be degenerate. 
An attempt to solve this problem has been carried out in a recent paper [IB], where 

an exactly separable version of the Bohr Hamiltonian, called ES-D, which uses a potential 
of the form u{p) + u{-f)/p'^ [39j, with a Davidson potential /?2 + p^/f^^ ^^^ [^ the 



place of u{P), and a steep harmonic oscillator centered at 7 = 0° as ^(7), is developed. All 
bands (e.g., ground, j3 and 7) in this model are treated on equal footing [50], depending 
on two parameters, the Davidson parameter (3q and the stiffness c of the 7-potential. The 
model is found [16] to be applicable only to well deformed nuclei (with i?4/2 > 3.0) due to 
the /3^ denominator in the ^(7) term. Nevertheless, it reproduces very well the bandheads 
and energy spacings within bands of almost all rare earth and actinide nuclei, with R4/2 > 
3.0, for which available data exists, as well as most of the B{E2) transition rates. The 
most glaring discrepancy concerns B{E2) values for the /3 band to ground band which are 
typically overpredicted by an order of magnitude. The two exceptions where ES-D does not 
provide a good description of energy spectra are ^^^Sm and ^^^Gd, which have previously 
been shown [221 EI] to be well reproduced with the infinite square well potential of the 
critical point symmetry X(5). Furthermore, the ES-D model provides insights regarding 
the recently found correlation [52] between the 7 stiffness and the 7-bandhead energy. 

Concerning the long standing problem of producing a level scheme with IBA SU(3) 
degeneracies within the framework of the Bohr Hamiltonian, in the framework of the ES-D 
model the spacings within the ground and (3 bands are identical, because of the oscillator 
term in the u{l3) potential. It is therefore enough to examine the conditions under which 
the 2"*", 4"*", G"*", ...levels of the /3 and 7 bands are degenerate. One is then led [16] to 
minimize the rms deviation between the even levels of the /3i and 71 bands for fixed value 
of the normalized /?i-bandhead, R0/2 = E{0~^)/E{2f). Numerical results indicate that a 
reasonable degree of degeneracy is obtained for L^ax = 10 and -R0/2 > 15, which is of 
physical interest, since experimental -R0/2 values extend up to 27 . In the case of ^^^Th, 
in which R0/2 = 14.8, one can see that the above mentioned experimental rms deviation 
between the even levels of the Pi and 71 bands is <Jjl^,y{L^ax = 10) = 1.142, while the 
corresponding theoretical quantity is cr^^^(Linax = 10) = 0.593 [IH]. Therefore, although 
the overall fit is quite good, the degree of degeneracy obtained from theory is less than 
the one indicated by experiment. One could conclude that the ES-D model does contain 
parameter pairs which correspond to approximate degeneracy of the low lying even levels of 
the (3i and 71 bands, while at the same time the spacings within the /3i band are identical 
to the spacings within the ground band, but the problem of reproducing a SU(3) spectrum 
from the Bohr Hamiltonian remains conceptually open. 

5 Conclusion 

The field of shape/phase transitions and critical point symmetries is rapidly expanding. 
Many additional references can be found in the review articles [T2|, [T3|, [53], as well as in 
Ref. [5i 
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